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Abstract 



We analyze the electron and neutron electric dipole moments induced by 
R-parity violating interactions in supersymmetric models. It is pointed out 
that dominant contributions can come from one-loop diagrams involving both 
the bilinear and trilinear R-parity odd couplings, leading to somewhat severe 
O I constraints on the products of those couplings. 

o 

O ■ I- INTRODUCTION 

o 

\ Fermion dipole moments provide severe constraints on possible new physics beyond the 

Standard Model. The foremost example is the electric dipole moment (EDM) of the electron 
or neutron whose current experimental bound is given, respectively, by |jl]J^ 



\de\ < 4.3 X 10~2^ecm, 

: < 6.3 X lO^^^^ecm. (1) 



In the Minimal Supersymmetric Standard Model (MSSM), apart from the CKM phase, CP 
violation can arise from additional complex parameters describing the soft supersymmetry 
(SUSY) breaking terms in the effective Lagrangian. Such complex soft parameters generate a 
light fermion EDM through one- loop diagrams involving sparticles (squarks/sleptons and/or 
gauginos). The current bounds on the neutron and electron EDM constrain the phases of 
soft parameters to be less than (9(10~^), or the sfermion masses to be larger than 0{1) TeV 
0] unless there are cancellations among various contributions . 

When R-parity conservation is not assumed, the MSSM allows additional lepton number 
(L) or baryon number (B) violating interactions. In this case, one usually invokes other type 
of (either unbroken or spontaneously broken) symmetry which would ensure that R-parity 
violating couplings are small enough to avoid too fast proton decay and too large neutrino 
masses 0. Imposing a symmetry which would forbid (or highly suppress) B-violating cou- 
plings while allowing non-zero L- violating couplings is an attractive possibility since then an 
interesting pattern of neutrino masses and mixing can be generated H. Furthermore these 
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couplings can provide a new source of CP violation, and thus the EDMs of light fermions 

m- 

The quark and lepton EDMs arising from the R-parity violating Yukawa couplings have 
been studied recently in Refs. To construct effective dipole moment operator, it is 

required to have not only the chirality-flip but also the insertion of even number of R-parity 
violating couplings. In other words, the resulting fermion EDMs involve some powers of AA 
and/or A A* as well as some fermion masses or chirality- flipping Yukawa couplings |^ where 
A denotes a generic R-parity violating Yukawa coupling. At one-loop level, EDM diagrams 
involving AA exist. However such diagrams are suppressed by additional Majorana neutrino 
(or sneutrino) mass which would compensates the L = 2 of AA in the EDM coefficients, and 
thus are negligibly small. It has been verified that important contributions can arise from 
two-loop diagrams involving the insertion of quark or lepton mass inside the loop and also 
the products of R-parity violating couplings such as AA* or A^A*^ 

In this paper, we wish to point out that sizable EDMs can arise at one-loop level in the 
presence of both the bilinear and trilinear R-parity violating interactions. Let us recall that 
generic R-parity violating superpotential includes the bilinear Higgsino-lepton mixing term 

W 3 fi,L,H2, (2) 

where Lj {i = 1, 2, 3) are the lepton SU{2) doublets, and H2 is the Higgs doublet with 
hypercharge Y = 1/2. There exist also R-parity violating scalar bilinears in the soft SUSY 
breaking potential: 

V3B,L,H2 + mljj^L,Hl (3) 

where Hi is the Higgs doublet with Y = —1/2. These bilinear scalar terms induce nonzero 
vacuum expectation values (VEVs) of sneutrino fields: 

(z>,)(xi?r(ff°)+mi*^^(i7?)- (4) 

The presence of these bilinear couplings yields a mass mixings between the R-parity even 
and odd particles, which are described by the sneutrino VEVs (z/j) and the Higgsino-lepton 
mixing parameter /Xj. Since they are R-parity odd and also fiips the fermion chirality without 
involving additional particles in the interaction vertex, these mass mixings are more efficient 
than the R-parity violating Yukawa coupling A for generating fermion EDMs. We will see 
that there exist one-loop diagrams for fermion EDMs which are proportional to the products 
of bilinear and trilinear R-parity odd couplings, e.g. (i/j)A or /i*A. Such diagrams are not 
suppressed by additional powers of small fermion masses or Yukawa couplings, and thus can 
provide strong limits on these product of R-parity odd couplings. 

The organization of this paper is as follows. In Section II, we define the basis on which our 
analysis is performed and summarize how to diagonalize the mixings between the R-parity 
even and odd particles which are induced by the bilinear R-parity violation. In Section HI, 
we classify and compute the main contributions to the quark and lepton EDMs involving 
both the bilinear and trilinear R-parity violating couplings. The calculational details are 
presented in the Appendix, and the bounds on couplings from the electron and neutron 
EDMs are summarized in Tables 1 and 2 under the assumption that those couplings take 
generic complex values. We conclude in Section IV. 
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II. PARTICLE MIXINGS INDUCED BY R PARITY VIOLATION 



For comparison with experiments, it is convenient to work in the quark and lepton 
mass eigenbasis. In this prescription, we leave the neutrinos in the charged lepton mass 
eigenbasis since the mechanism of neutrino mass generation is quite different from that of 
other fermions. R-parity violation under consideration can be a major source of nonzero 
neutrino masses, and then the charged lepton mass eigenbasis is particularly useful for 
analyzing neutrino mass matrix. In the quark and charged lepton mass eigenbasis, the full 
superpotential of the MSSM fields including generic R-parity violating couplings is given by 



where Li and Qi are the lepton and quark SU{2) doublets, and f/f, are the SU{2) 
singlet lepton and anti-quark superfields. Here Vij is the CKM matrix of quark fields, and 
the R-parity violating lepton-Higgsino mass mixing is described by dimensionless parameter 
ej in the unit of the conventional Higgsino mass ^. From the above superpotential, one 
obtains the following lepton number violating Yukawa vertices in the Lagrangian; 



+ Kjk (j^idkRdjL + djLdkR^iL + d^^^vf^djij 

- Kik^ij {endkRUji + UjidkR^iL + dlRemUji) , 

where the superscript C denotes the charge conjugation. In order to ensure the longevity 
of proton, the products A'A" have to be highly suppressed [|10|- Bilinear R-parity violating 



terms can also contribute to proton decay through the products {vijX" or {eijj,)*X" |n]]. In 
this paper, we simply assume that A" is small enough to ensure the proton stability, while the 
L-violating couplings can take arbitrary values as long as they satisfy other phenomenological 
constraints. 

So far we have defined the quark and lepton mass eigenbasis without including the 
effects of R-parity violation. In the presence of bilinear R-parity violation, there appear 
mixings between the ordinary fermions (Higgses) and the gauginos/Higgsinos (sleptons). As 
a result, further diagonalization is required to define the true mass eigenstates. In the below, 
we present the mixing matrices between the R-parity even and odd particles relevant for 
our calculation of EDMs. For this, we generalize the results in Ref. |jl2| to take into account 



the complexity of the bilinear parameters in Eqs. (^,|^,^), which leads to a mixing between 
the CP even and CP odd neutral Higgses and sneutrinos. The process of diagonalization 
of slepton-Higgs mixing involves a rotation by which the Goldstone modes in the Higgs- 
slepton sector are decoupled |T^. For the sake of simplicity, we go to the new basis in 



which the lepton and Y = —1/2 Higgs superfields are redefined so that sneutrinos have 



vanishing VEVs |13|. When the sneutrino VEVs in the original basis are small enough, i.e. 
Qi = {u*) / {H^) <^ 1 in Eq. (0), the Lagrangian in the new basis can be obtained by the 
field redefinition 
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Li = Li + a* Hi , Hi = Hi- aiLi . (6) 
Then, the R-parity odd parameters in the new basis are given by 
Bi = Bi- QiB , ml^H, = ^L^H, + aii^l, - ^h,) 

h = — O-i 1 Kjk = Kjk — O'ih'^jSjk , ^{jk = Kjk ~ O'ih'^j^jk (7) 

where /i^ and are the Yukawa couphngs of the charged leptons and down-type quarks, 
respectively. Note that we have the relation I3itp+rh\,^^ = where = tan P = {H2) / {H^), 
which assures that sneutrino VEVs vanish in the new basis. In the following, we will drop 
hats to denote quantities in this new basis of vanishing sneutrino VEVs. 

Neutrino-neutralino mixing: Let us first consider the well-known mixing among 
neutral fermions for the completeness. In the basis of vanishing sneutrino VEVs, the 
mass mixing between the neutrinos (z/j) and the neutralinos Nj = [B, W3, H^, H2) are 
described by the following 3x4 matrix: 

(MD),, = (0,0,0,-e,/x) . (8) 

Combined with the usual 4x4 neutralino mass matrix Mn of Nj = {B,W3, H^, H2), the 
above 3x4 mixing matrix Md forms a 7 x 7 mass matrix of the neutrinos and neutralinos. 
At leading order in the expansion in powers of small Md, the 7x7 neutral fermion mass 
matrix can be diagonalized by the rotation: 



" =U(e«r 1 jlo ivj 

where 6^ = M^iM^^ U"" is the 3 x 3 neutrino mixing matrix, and finally is the usual 
4x4 diagonalization matrix of Mn in the absence of R-parity violation Here the 

mixing matrix 0^ is given by [O 



= e,cf + e,5,3 (10) 



with 



^ Mz M2 Ml Mz Mz\ 

where Fn = MiM2/M^-Mls2p/ ^i, = c^Mi-Fs^Ma for the SU{2) x ?7(1) gaugino mass 
parameters Mi, M2, sw = sin 9w with 9w being the weak mixing angle, and S/3 = sin/3, etc. 
In the above, for the sake of simplicity, we assume that Mi, M2, and fi are approximately 
real, which may be necessary also to avoid a too large neutron EDM. 

Charged lepton-chargino mixing: R-parity violation leads also a mixing between the 
charged leptons and the charged gaugino /Higgsino. Following the steps similar to the case 
of neutral fermions, one finds that the 5x5 diagonalization matrices U^,U^ which would 
bring (e,, W~ , Hi ) and (e^, W~^, H2 ) to the mass eigenstates (ef , xf , xf) are given by 

,L.B ( 1 <S>''"\(\ 
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where 



©fa = e.CC^ + ^i5a2 , ©fa = TT^QCC/? 



with 



(13) 



- U2 



VI/ 



^ ^ _ / ^ Mwifi-M^tp) (M| + 2M^c 

Here Fq = M2 + My^si^/ and mf denote the charged lepton mass. Note that 0^ is 
suppressed by the factor my Fc compared to 6^. In Eq. (|12D , the 2x2 matrices L, R are 
the usual chargino mass diagonahzation matrices in the absence of R-parity violation ||14|| , 
and the index a = 1, 2 labels the chargino fields before the diagonahzation by L, R. 

Charged slepton-charged Higgs mixing: In the basis with vanishing sneutrinos 
VEVs, no mass mixing between the slepton and Higgses arises from the D-term potential. 
Then the mixing comes from the following part of the scalar potential: 



V 3 Bi{UH2 - tpUHl) + eiht/jH^HlE^ + h.c, 



(14) 



where we have used the relation mj^.fj^ = —Bitp for vanishing sneutrino VEVs. Decoupling 
the Goldstone mode by setting iff = S/^if" and H2 = cpH~ , we find that the mixing mass 
terms for the charged slepton and Higgs fields are given by 



V 3 -{Bi/cp)eiLH+ - {m1pi/cp)eiRH^ + /i.e., 



(15) 



where en and em are the left-handed and right-handed sleptons, respectively. Applying 
again the seesaw diagonahzation as in the case of Eqs. (P,0), one finds the 8x7 rotation 
matrix which brings the charged scalar fields (e^x, Cj/j, iff, iff) to the mass eigenstates 



U 



c 



( 1 


-ri\*Sfs 




1 



ri\\ 



(16) 



where 



V2 



for M|-, M]^- and M|,^ being the left-left, right-right and left-right mass terms of slepton fields 
in the case of R-parity conservation, respectively. Here we have neglected the diagonahzation 
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of the slepton mass matrix which is assumed to be an approximate unit matrix for reasonable 
range of tan (3. 

Sneutrino-neutral Higgs mixing: If we allow fn^.^^ and/or Bi to be complex- valued, 
these R parity violating parameters lead to a mixing not only between the R-parity even and 
odd states, but also between the CP even and odd states. Then the usual decomposition 
of slepton/Higgs sector into CP-even and CP-odd parts [|l^ cannot be applied any more. 
Again decoupling the Goldstone mode by setting \m.{Hf) = SjsA^ and lm{H2) = cpA^, we 
find the mixing mass terms for (Re(i>j), Im(z/j)) and {R.e{Hi),Iie{H2), A^) are given by 



Re{Bi) 
-MBi 



-lm{Bi)/cp\ 
-Re{B,)/cJ 



(17) 



Combining the above mixing mass matrix with the R-parity conserving slepton and Higgs 
mass matrix, we get the diagonalization matrix given as follows: 







(18) 



where the index i labels each slepton generation and Ua is a block-diagonal 3x3 matrix 
in which the upper-left block is a 2 x 2 matrix diagonalizing the CP-even Higgs bosons in 
the absence of R-parity violation |T^ and the lower-right block is identity. Here we keep 
using the standard notation for the Higgs-like scalar fields although they are not exact CP 
eigenstates. Again, the slepton diagonalization part is neg lected. In Eq. ([T|), (for each 
slepton flavor z) is a 2 x 3 matrix whose components are given by 



^11 



^12 



[m\ + M|c2/^)t^Re(E,) 
Fs 

{m\ - M|c2/5)Re(i?.) 



s _ Im(A 



'^21" 



e 



22" 



m 



2 ' 



e 



{ml + Mlc2p)tp\m{Bi 
Fs 

(m? - M|c2^)Im(i?,) 



Re{Bi 



(19) 



23" 



where F< 



s 



-m 



■m'jjo)+'m'loTn'j^o and m?. = m|^. + M|c2/3/2. As we will see in the 
Appendix, the EDMs arising from slepton-Higgs mixing turns out to be proportional to the 
charged lepton or down-type Yukawa couplings [see the terms containing in Eqs. ( |A2| , |A3D ] 
and thus are subleading. 

An important consequence of neutrino-neutralino mixing is the generation of nonzero 
neutrino mass Diagonalizing the 7x7 neutral fermion mass matrix, one finds the 
neutrino mass matrix at tree-level: 



m. 



tree 



Ml 



(20) 



The existing atmospheric and solar neutrino data suggest that the largest neutrino mass 
eigenvalue is of the order of 0.1 eV [|T6|. Also the neutrino masses induced by R-parity 
violation are generically hierarchical and dominated by the above tree level mass matrix. 
We then have 
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where niy^ denotes the largest mass eigenvalue of (PO). If R-parity violation is the dominant 



source of neutrino masses, the large mixing of atmospheric neutrinos ||T6l would imply e ~ 
1^2! ~ l^sl, and one can have |ei| ^ e or |ei| ^ 0.04e depending on the choice of the large or 
small mixing angle solution to the solar neutrino problem |TT 



III. FERMION EDM FROM R-PARITY VIOLATION 

The fermion EDM operator is given by 

'-d^^a^^'^^^iPF,, (22) 

for a Dirac fermion and the electromagnetic field strength F^^. Obviously this EDM 
operator flips the fermion chirality and also breaks CP invariance. It has been argued that 
there is no one-loop fermion EDM from R-parity violating trilinear Yukawa couplings in 
the limit of massless neutrino J^,^- To generate R-parity even EDM, one needs to insert 
even number of R-parity odd couplings. Then it is not possible to have one-loop EDM 
diagrams involving the two insertions of R-parity violating Yukawa couplings in a manner 
compatible with both the requirement of chirality-fiip and the chiral structure of Yukawa 
vertices. However in the presence of bilinear R-parity violation which would yield a mixing 
between the R-parity even and odd particles, a combination of R-parity violating bilinear 
and trilinear couplings can induce a fermion EDM at one-loop level. An important feature 
of these one-loop diagrams is that the second trilinear interaction vertex is provided by 
gauge interactions (or the large top quark Yukawa coupling), so there is no small fermion 
mass or Yukawa couplings involved other than the minimal product of R-parity odd bilinear 
and trilinear couplings. FIGs. 1-3 show such diagrams for the EDMs of charged lepton, 
down-type quarks, and up-type quarks, respectively. One can also have one-loop diagrams 
involving the two insertions of R-parity odd bilinear couplings. However such diagrams are 
suppressed by small neutrino mass, and turn out to be negligibly small. 

Let us first consider the leading contributions to the electron EDM in more detail. The 
diagrams of FIG. 1 without a flavor change of sleptons inside the loop, i.e. i = k = a = \, 
give rise to the electron EDM : 

+V2tw Im «„Wj^„Ai,i) -^a(m,o; mA , (23) 

where Gf, Gg are the loop functions defined in the Appendix and are the diag- 

onalization matrices discussed in the previous section. Here Xni Xn ^'^^ fermions inside 
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the loops of FIG. 1 which can be either R-parity even, i.e. charged lepton or neutrino, or 
R-parity odd, i.e. chargino or neutrahno. It turns out that the terms with R-parity even 
Xn is negligible, so we will keep only the terms with R-parity odd Xn- In this case, l^wxn^ 
Uwxn^ ^Bxn correspond to the R-parity even elements of the usual chargino or neutralino 
diagonalization matrix, while W^^^ or Wj^^^ contains the R-parity odd element 6^ or 
which connects the neutrino with neutralino or the charged lepton with chargino. 

Applying the experimental bound on the electron EDM to Eq. (P5D, we can obtain the 
bounds on the products of the trilinear and bilinear R-parity violating couplings. For a 
rough estimate, it is useful to recall that the mass mixing between the charged wino 
and the charged lepton ej (denoted by the cross in the fermion line of FIG. la) is given by 
MwejC/3, implying 

Then applying the experimental bound (|I]) to the first term of Eq. (pSf), we get 

Im (A,,.-c,) . Im (a,,.-^) < 1.4 x lO^^^^ . (24) 

Here the bound on the product Aijie*Mvi//m^± can be understood by that the mixing mass 
of H2 and the charged lepton is given by e^/i and the Wino-Higgsino mixing is of the 
order Mw/m-^ with Mw < m^. For numerical estimates, we took a typical value of the 
loop functions, G f{mf;ma) ~ 1/3 and Gs{mf] nia) ~ 1/6. Note that the values of these 
loop functions are rather insenstive to mf/rria as can be seen in FIG. 4. The other terms 
in Eq. ( |23D are expected to yield a bit weaker bound by the factor of tw If the dominant 
source of neutrino masses were the bilinear R-parity violation, we have ec/j = (F^mj^Y^'^ / Mz 
as discussed in the previous section. This enables us to translate the bounds of Eq. ( |2^ 
into a bound on the R-parity violating Yukawa coupling: 

where we assumed that the phases are of order unity and F/v ~ m^± . 

If we allow a generation mixing in the squark or slepton sector, we have more contribu- 
tions to the fermion dipole moments. To deal with this generation mixing, let us use the 



so-called mass insertion method [|l^]. In this approach, the amount of flavor violation is 
parameterized by the quantity, 

Sif' ^ , (26) 

1717 771? 

where Afj^^ is the off-diagonal element and is the diagonal element of the sfermion 
mass-squared matrix in the quark and lepton mass eigenbasis, f = u,d,l are the indices for 
the up-type squarks, down-type squarks, sleptons, and A, B = L, R denote the chirality of 
these sfermions. Including the flavor mixing of sfermions, we have for instance the following 
contribution to the electron EDM: 
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eg 



(27) 



which reproduces the first term in Eq. (|23|) when 1 = a. The definition of the loop function 
Hf and also the other contributions can be found in the Appendix. Again applying the 
experimental bound to these contributions to the electron EDM, we get 



) , Im (4''A.,ie*^) < 1.4 x 10 



(lOOGeV)^ 



Im (S^ia^Xijo,e*C/3) , Im 6^{^^Xijae", 



w 



< 6.4 X 10"® ^ ' 



(100GeV)2 



(28) 



where the first bound comes from the diagram of FIG. la or lb exchanging u or e and the 
second one from FIG. Ic exchanging e'^. Most stringent bounds on d^'j with AB = LL, RR 
come from the /j Ij'y decays. One then finds Si2^ ^ 10"^(m^-/100 GeV)^, however there 



are essentially no bounds on 6i^^ and 62^^ [ p!^ . 

The most important contributions to the down quark EDM come from the gauge vertex 
diagrams [FIG. 2a, 2b], and also the top Yukawa vertex diagrams [FIG. 2c], from which we 
find 



:1,AB 



eg 

167r2 
1 



1 

+ — %Im (w^^^Wj^J^Aj'ii) ^a(m^o;mj.J 



eht 
167r' 



■ImfwgVW^-'^aiVjU' 



(29) 



More general expression including the squark mixing effects is given in the Appendix. From 
the valence quark contribution to the neutron EDM 



dn = -{4:dd - du^ , 

taking into account the renormalization group effect from the weak scale to the hadronic 
scale with r*^ ~ 1.5 [|T^], the current experimental bound (|I|) gives 



Im(KiV^lA;,i6*c^,) , Im [v^^X'.n^^ 



w 



< 1.0 X 10"^ ^ ^ 



(lOOGeV)^ 



(30) 



Again if the bilinear R-parity violation gives my ^ 0.1 eV, the above bound would give 
■^jii ~ O-l- From the top quark exchange diagram [FIG. 2c] giving the last contribution of 
Eq. (P^D, one finds 
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as can be read off from Eq. (0). The contribution involving is neglected as it is propor- 
tional to the lepton mass m*. Then one gets 



< 2.5 X 10" 



0.01 



rrir 



Vsi 100 GeV 



(31) 



If we allow the squark generation mixing, VaiVil^X'ji^ can be replaced by A^q,]^5^{^^ in Eq. (^U|). 
More general expressions with both the CKM and squark mixings are given in the Appendix. 

Contrary to the down quark EDM, the up quark EDM in our scheme turns out to be 
very small as it is further suppresed by small quark masses. For instance, the diagram in 
FIG. 3 gives rise to 



eh 



mi 



3 



m - ; m 



+ G, 



S'^d ) 



(32) 



which is the limiting case of a = /3 in Eq. (^J). One can see from Eq. (|32D that the up 
quark EDM includes additional suppression factor m^m^/m?. 

The nucleon EDM gets contributions not only from the quark EDMs, but also from the 
chromoelectric dipole moments (CDMs) of light quarks [^. Since the diagrams for the 
CDMs of the up and down quarks have the same structure as those for the EDMs under the 
replacement of the external photon by gluons, one does not find a new bound on the R-parity 
violating couplings. On the other hand, there is also a contribution from the strange quark 
CDM whose diagrams involve a different set of R-parity violating couplings. The strange 
quark CDM (rf^) contribution to the neutron EDM is given by dn ~ 0.027(igr^ ||21|] with 
r'^ ~ 3.3 coming from the renormalization group effect, which leads to 



Im (K.2V;lA^;2e*C/3) , Im \42VJ^A^.j2e* 



M, 



w 



< 1.7 X 10 



_g TTt-y/VTiiq^ 

(100GeV)2 



(33) 



Collecting all the contributions to the electron and neutron EDMs, we summarize the re- 
sulting bounds on R-parity violating couplings in the Tables. 



IV. CONCLUSION 

We have analyzed the electron and neutron EDMs induced by R-parity violating inter- 
actions, taking into account both the bilinear and trilinear R-parity violating terms. It is 
pointed out that the fermion EDMs can arise at one-loop level as a combined effect of the 
bilinear and trilinear terms. This has to be contrasted to the previous analyses showing 
that nonzero EDMs arise at two-loop level considering only the trilinear R-parity violating 
vertices. We have computed the complete one-loop diagrams (except those proportional 
to tiny neutrino masses) involving both the bilinear and trilinear parameters as presented 
in the Appendix. Among them, the leading contribution comes from the diagrams with a 
gauge/top- Yukawa vertex and a R-parity violating trilinear vertex for which the chirality- 
flip is provided by (bilinear) mass mixing between R-parity even and odd particles. These 
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one-loop diagrams lead to somewhat severe constraints on the products of the R-parity vi- 
olating bilinear and trilinear couplings as summarized in Table I. If the bilinear R-parity 
violation is a dominant source for the neutrino mass rrii, ~ 0.1 eV as was suggested by the 
Super-Kamionkande atmospheric neutrino data, the typical upper bounds of O(10~^) can 
be put on some Xijk couplings involving the first generation index as shown in Table 11. 

Acknowledgement This work is supported by Seoam Foundation (K.C.) and by grant 
No. 1999-2-111-002-5 from the interdisciplinary Research program of the KOSEF and BK21 
project of the Ministry of Education (K.C. and K.H). 



11 



REFERENCES 



[1] K. Abdullah et al, Pliys. Rev. Lett. 65, 2340 (1990); E. Commins et al, Pliys. Rev. A 50 
(1994) 2960; B.E. Sauer, J. Wang and E.A. Hinds, Phys. Rev. Lett. 74, 1554 (1995); J. 
Chem. Phys 105 (1996) 7412. 

[2] P.G. Harris et al, Phys. Rev. Lett. 82, 904 (1999). 

[3] J. EUis, S. Ferrara and D.V. Nanopoulos, Phys. Lett. B114, 231 (1982); W. Buchmuller 
and D. Wyler, Phys. Lett. B121, 321 (1983); J. Polchinski and M.B. Wise, Phys. 
Lett. B125, 393 (1983); J-M. Gerard, et al, Nucl. Phys. B253, 93 (1985); M. Dugan, 
B. Greenstein and L. Hall, Nucl. Phys. B255, 413 (1985); A. Sanda, Phys. Rev. D32, 
2992 (1985). 

[4] T. Ibrahim and P. Nath, Phys. Lett. B418, 98 (1998); Phys. Rev. D57, 478 (1998); 
T. Falk and K.A. Olive, Phys. Lett. B439, 71 (1998); M. Brhlik, G.J. Good and 
G.L. Kane, Phys. Rev. D59, 115004 (1999); S. Pokorski, J. Rosiek and C.A. Savoy, 
Nucl. Phys. B570, 81 (2000). 

[5] V. Ben-Hamo and Y. Nir, Phys. Lett. B339, 77 (1994); T. Banks, et al, Phys. Rev. 
D52, 5319 (1995); E. J. Chun and A. Lukas, Phys. Lett. B387, 99 (1996); K. Choi, E. 
J. Chun and H. Kim, Phys. Rev. D55, 7010 (1997); K. Choi, E. J. Chun and H. Kim, 
Phys. Lett. B394, 89 (1997); G. Eyal and Y. Nir, JHEP 9906, 025 (1999); Q. Shafi and 
Z. Tavartkiladze, |hep-ph/9909"238i 

[6] L. Hall and M. Suzuki, Nucl. Phys. B231, 419 (1984); F. M. Borzumati, et al, Phys. 
Lett. B384 (1996) 123; B. de Carlos and P. L. White, Phys. Rev. D54 (1996) 3427; 
A. Yu. Smirnov and F. Vissani, Nucl. Phys. B460 (1996) 37; R. Hempfling, Nucl. Phys 
B478 (1996) 3; H. P. Miles and N. Polonsky, Nucl. Phys. B484 (1997) 33; E. Nardi, 
Phys. Rev. D55 (1997) 5772; E. J. Chun, et al, Nucl. Phys. B544, 89 (1999); A. S. 
Joshipura and S. K. Vempati, Phys. Rev. D60, 095009 (1999); J. Ferrandis, Phys. Rev. 
D60, 095012 (1999); M. Bisset, et al, |hep-ph/981149q ; S. Rakshit, G. Bhattacharyya, 
Phys. Rev. D59, 091701 (1999); R. Adhikari, G. Omanovic, Phys. Rev. D59, 073003 
(1999); K. Choi, E. J. Chun and K. Hwang, Phys. Rev. D60, 031301 (1999); D. E. 
Kaplan and A. Nelson, |hep-ph/990T25l A. Abada and M. Losada, |hep-ph/990835^ ; J. 
L. Chkareuh, et al, |hep-ph/9908451| ; J. C. Romao, et al, |hep-ph/990749^ ; O. Haug, et 
al, |hep-ph/99093T8| . 



[7] R.M. Godbole, S. Pakvasa, S.D. Rindani and X. Tata, [hep-ph/9912315 



[8] S.A. Abel, A. Dedes and H.K. Dreiner, |hep-ph/9912429 



[9] S.M. Barr, E.M. Friere and A. Zee, Phys. Rev. Lett. 65, 2626 (1990). 
[10] L Hinchliffe and T. Kaeding, Phys. Rev. D47, 279 (1993); A.Y. Simrnov and F. Vissani, 

Phys. Lett. B380, 317 (1996). 
[11] B. Mukhopadhyaya, et al, Phys. Lett. B443, 191 (1998); E.J. Chun and J.S. Lee, Phys. 

Rev. D60, 075006 (1999). 
[12] E.J. Chun and S.K. Kang, Phys. Rev. D61, 078012 (2000), ||hep-ph/ 9909429 . 



[13] Y. Grossman and H.E. Haber, Phys. Rev. D59, 093008 (1999). 
[14] For a review, see H. Haber and G. Kane, Phys. Rep. 117, 75 (1985). 
[15] A.S. Joshipura and M. Nowakowski, Phys. Rev. D51, 2421 (1995); M. Nowakowski and 
A. Pilaftsis, Nucl. Phys. B461, 19 (1996). 



12 



[16] The Super-Kamiokande Collaboration, Y. Fukuda, et ai, Phys. Rev. Lett. 81, 1562 
(1998). 

[17] For recent analyses, see, J.N. Balicall, P.I. Krastev and A.Y. Smirnov, Phys. Rev. D58, 
096016 (1998). 

[18] F. Gabbiani, et al, Nucl. Phys. B477, 321 (1996); M. Misiak, S. Pokorski and J. 
Rosiek, in "Heavy Flavours 11", ec? A.J. Buras and M. Lindner (World Scientific, 1998), 
||hep-ph/9703442|| . 

[19] E. Braaten, C.S. Li and T.C. Yuan, Phys. Rev. Lett. 64, 1709 (1990); R. Arnowitt, J. 

Lopez and D.V. Nanopoulos, Phys. Rev. D42, 2423 (1990). 
[20] X.G. He, B. McKeller and S. Pakvasa, Int. J. Mod. Phys. A4, 5011 (1989); Erratum, 

ibid. A6,1063 (1991) 
[21] X.G. He, B. McKeller and S. Pakvasa, Phys. Lett. B254, 231 (1991). 



13 



TABLES 



coupling 


upper bound 


Xijie*cp 

Ol3 ^231620^ 

rl.LL \ d; 


1.4 X IQ-^Xs 
1.4 X 10-^Xs 

U.D4 X lU Xg 

0.64 X lO-^a;^ 


^l"^«2A;i2S*C/3 


1.0 X IQ-'^Xs 
1.0 X lO^^Xs 
1.7 X lO^^Xs 
1 n-6 0.01 \ 

/.D X iU lOOGeV j 


TABLE I. Bounds on the imaginart part of the combinations of trihnear and bilinear couplings. 
Note that e^c^ can be replaced by (*jM\\i jm^. Here Xg = (loo^eV^' 


coupling 


upper bound 


Al21, Ai3i 

^Yi A231 

A,LL , 
Ol3 ^^123 

Ol3 ^133 


O.OlSx^x,^ 

Qmzxsx,, 

Q.Q^^XgXi, 
Q.QhQxgXi, 


^13 '^jai 


O.lOx^jXjy 
0.45X3X1^ 
0.10X5X1/ 



TABLE II. Bounds on the magnitude of single trilinear coupling under the assumption that 
neutrino masses come dominantly from the R-parity violating bilinear terms. As a reference, we 
take m^, - 0.1 eV and then x„ = {mOSv)'^ imSv)~^ ■ 
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APPENDIX A: FULL ONE-LOOP CONTRIBUTIONS TO THE FERMION 
DIPOLE MOMENTS FROM R-PARITY VIOLATION 



In this Appendix, we provide complete expressions for the one-loop fermion dipole mo- 
ments induced by R-parity violating couplings. For this, let us first define some loop func- 
tions. For the basic loop functions given by 



Gf{t) 
Gsit) 



t 



2(1-^)2 {l-tf 
t + 1 tint 



the modified loop functions which are propotional to the chirality-fiipping fermion masses 
are defined as 



Gf{mf,m^) = —Gj 



rrin 



GsimKrUa) = —Gs 



m 



f 



ml 



These functions behave well in the limit m/ — > and mildly depend on mf/rria as shown in 
FIG. 4. When the mass insertion approximation is used, another useful set of (dimensionful) 
loop functions are 



Hf{mf;ma,'mp) 
Hs{mf;ma,mp) 



ml — m"^ 
mfmamp 



-Gi 



mf, 



mt 



m 



— kGs 
ml 



m. 



mf, 



m. 



mt 



m} 

2 "-^^ 



^ mY 
'mj^ 



Any contribution to the electric dipole moment between the two fermions fi and fk can be 
expressed in the following form: 



d 



ki- 



(Al) 



where / stands for the external quarks or charged leptons, i, k are their generation indices, 

AB denote the chirality of intermediate sfermions, i.e. AB = LL, LR, RR, and a and f5 
are the indices for those sfermions. Here we do not seperately treat the diagrams without 
sfermion flavor mixing. The expressions for such diagrams can be easily obtained from those 
for the diagrams with sfermion flavor mixing by changing Sj^p^, Hf, Hg as follows: 



Hf{mf;ma,mf3) 
Hs{mf;ma,mf)) 



— Gf{mf;ma) , 



m, 



— Gs{mf;ma) ■ 
mry 
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by 



The part {S ■A)ik for the charged lepton dipole moments between li and Ik are then given 



(A2) 



I "''4>i ) 



A 



jik 



+ htl^H^xn^^Jai^^jakHs{m^o^;me„,m~c.) 

+ K^HiXn^fxn^ak^^jiocHs{m^jm~c^, m~e^) 

and for the dipole moments between the down type-quarks di and dk, 



(A3) 



9 



\ ^4>i ) 



X' 



jik 



+ \hPlxMfxJT>^3Mrn^o^-,m^^,rn^.) 



and finally for the up type quarks itj and Uk, 

{S . A)l = - /7^Z^g,-Wf,-l^«KJ,AL,^ {-lGf{md,;m,-) + G,(m,,; m,-)} (A4) 

Here 0f denotes the charged and neutral scalar mass-eigenstates, respectively. 
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FIGURES 



9 W~ 



\ 

\ 



-X- 



ei g',g B,W3 



^ajk ^k 



(b) 



/ \ 

/ \ 
/ \ 
/ \ 
/ \ 
/ I 

^ . X ^ ' — - 

ei -\ja ^3 B g' 



FIG. 1. Most important diagrams for the charged lepton EDM. The black blobs denote possible 
scalar mixing and the crosses are chirality-fliping mass insertions. 
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Ur 



M/3 



-X- 



(a) 



-X 



(]. X'.. Uj 



B g' d% 



(b) 



/ \ 
/ \ 
/ \ 
/ \ 

I 1 



-X 



d^ hVs. t a;,,\/,^ dl 



FIG. 2. Most important diagrams for the down-type quark EDM. There exist also diagrams 
similar to FIG. lb, and also a supersymmetric counterpart of (c). 



18 



-X- 



FIG. 3. Leading diagram for the up-type quark EDM. There is also a supersymmetric counter- 
part of this. 
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FIG. 4. Loop functions Gf{mf;ma) and GsimfjUia) as a function of t = m// 
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